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; Abstract. The source of Ultra High Energy Cosmic Rays (UHECR) is still an unresolved 

mystery. Up until recently, sources of Gamma Ray Bursts (GRBs) had been considered as a 
suitable source for UHECR. Within the fireball model, the UHECR produced at GRBs should 
be accompanied with a neutrino flux detectable at the neutrino telescope such as IceCube. 
q ■ Recently, IceCube has set an upper bound on the neutrino flux accompanied by GRBs about 

CN . 3.7 times below the prediction. We investigate whether this deficit can be explained by the 

oscillation of the active neutrinos to sterile neutrinos en route from the source to the detectors 
within the pseudo-Dirac scenario. We then discuss the implication of this scenario for diffuse 
^ . supernova relic neutrinos. 
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1 Introduction 

One of the challenges in cosmic ray physics is identifying the source of Ultra High Energy 
Cosmic Rays (UHECRs) with energy > 10 18 eV. Simple calculation of the energy budget 
and rate of astrophysical explosions shows that the transient Gamma Ray Bursts (GRBs) 
can accommodate the observed flux of UHECRs [1]. The GRBs are energetic explosions 
with luminosity ~ 10 51 erg/s at cosmological distances which have been observed by various 
satellites via their gamma ray emissions at a rate ~ 2/day. 

Within the fireball model of the GRB explosion (see [2] for a review on this model), 
the decay of the photo-pions produced by accelerated protons in the expanding shock wave 
generates a flux of accompanying neutrinos. Based on the observed flux of UHECRs, the 
expected flux of neutrinos is EldN u /dE u ~ 5 x 10~ 9 GeVcm _2 s _1 sr _1 in the energy range 
~ 100 TeV — lOPeV (the so-called Waxman-Bahcall flux [3]). On the other hand within 
the fireball model, the neutrino flux emitted from each GRB can be calculated from the 
observation of gamma ray energy spectrum. Normalizations based on the cosmic ray flux and 
the gamma ray flux both result in approximately the same prediction for the neutrino flux. 
The Waxman-Bahcall flux leads to ~ 10 events per year in a km 3 -scale neutrino telescope. 

Recently, the IceCube experiment published the result of the analysis of data taken 
during the construction of detector, that is IC40 and IC59 [4]. The combined analysis of 
the data does not show any neutrino signal correlated with the observed GRBs during the 
data-taking period. The limit of IceCube on the neutrino flux is ~ 3.7 times smaller than 
the prediction [4]. From this limit, the IceCube collaboration has concluded that the GRBs 
cannot account for the bulk of the UHECRs in the universe [4]. However, this conclusion can 
be questioned in three ways: (1) As shown in [5], within the fireball model, uncertainties are 
so large that the prediction for the neutrino flux can be lower than the standard Waxman- 
Bahcall prediction by up to a factor of ten. (2) In [6] , it is discussed that within the cannonball 
model, the sources of GRBs can account for the UHECR without violating the neutrino flux 
bound because the cannonball model does not predict a sizable neutrino flux. (3) Finally, 
it is possible that some beyond standard model effects reduce the flux after production (see 
e.g. [7]). In this letter, we entertain the third possibility. 

Although overwhelming evidence has been gathered showing that there are at least two 
massive neutrino mass eigenstates, we do not still know whether neutrinos are of Majorana 
type or of Dirac type. An interesting possibility is to have pseudo-Dirac neutrino mass 
scheme [8] , where the dominant contribution to the neutrino mass comes from the Dirac mass 
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term, mc, with a small correction from the Majorana mass term, tum such that mu <C ttid- 
Notice that the presence of the Dirac term requires right-handed sterile neutrinos. Within 
the pseudo-Dirac scenario, there is a small mass splitting between active and sterile neutrinos 
making the oscillation between active and sterile components in principle possible. For very 
small ttim , the oscillation length will be too large for the atmospheric and solar neutrinos 
to oscillate to sterile ones. However for neutrinos traveling over cosmological distances, the 
baseline will be large enough for active-sterile oscillation to take place [9]. In this paper, we 
investigate this possibility to relax the conflict between IceCube bound on the neutrinos from 
GRBs and the expected neutrino flux from GRBs as the origin of UHECR. 

The paper is organized as follows: In section 2, we briefly discuss the pseudo-Dirac 
scenario for neutrino masses. In section 3, we investigate the possible depletion of GRB 
neutrino flux due to the presence of tiny active-sterile mass-squared difference predicted in 
pseudo-Dirac scenario. Section 4 is devoted to the implications of pseudo-Dirac scenario for 
supernova neutrinos, discussing both single source and diffuse fluxes. The conclusions are 
summarized in section 5. 

2 Pseudo-Dirac Scenario 

In general, a neutral fermion such as neutrino can have two types of mass term: (1) Dirac 
mass term, mr>vv\ (2) Majorana mass term, ^m^zZ^Cf — i/^Cu*), where C is the charge 
conjugation operator. Dirac mass term requires the existence of a right-handed neutrino: 
mvvv = mnipRVL + ^lvr)- If there is only Dirac mass term, the neutrino mass is said 
to be purely of Dirac type and the lepton number is conserved. For the case of pure Dirac 
neutrino mass term, the mass eigenstates will be the superposition of Majorana states with 
maximal mixing and the mass eigenvalues will be equal. To illustrate this point, consider one 
generation pure Dirac mass term. In this case, the two Majorana states u + = (vl + 
and v~ = — i^|)/v2 will be mass eigenstates with masses equal to m + = mp and 
m~ = —rriD, respectively. The two states vr and vl are degenerate in mass, which taking 
into account that their masses are lighter than implies that vr has to be sterile. Notice 

that despite the maximal active-sterile mixing, there will be no oscillation between active and 
sterile neutrinos because the squares of the masses are exactly equal. 

Let us now suppose, through a lepton number violating process, the right-handed neu- 
trinos obtain a very small Majorana mass, such that tum ^ mD (which justifies the nomen- 
clature "pseudo-Dirac"). The smallness of tum can be justified within the 't Hooft criterion 
as in the limit mu — > 0, the lepton number is conserved. In this case the degeneracy between 
(m + ) 2 and (m, - ) 2 will be lifted and the active-sterile oscillation can take place. Generaliza- 
tion to three generations of neutrinos is straightforward. In this case the generic mass term 
consisting of both Dirac and Majorana mass terms, can be written as C mass = — ^^ C M^, 
where \& = {vlXi u L2, ^L3, v bx^ v< r2i v m) T ^ an d the 6x6 mass matrix is given by 



where mo and vtim are the 3x3 Dirac and right-handed Majorana mass matrices, respectively. 
It is convenient to go to the mass basis where mo is diagonal: 





mo = diag (m 1 ,m 2 ,m 3 ) 
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where mi = y mf + Am^ and = y m\ + Amf tm , with Am;: ol and Am^ tm , respectively, 
being the solar and atmospheric mass-squared differences obtained in neutrino oscillation 
phenomenology. In this basis, ium is a general symmetric matrix with mu <C mu (order 
relation applies to the nonzero elements of the matrices). The mass eigenstates are given by 



VLi + V 



a* 



7r ^ + EK^' + ^ 



(2.2) 



where the coefficients of- , j3, 



± 



mM/ m D *C 1 and the following orthogonality relations 
apply to the mass eigenstates: (z/^ |za~) = and (v^\vj) = (Vi\Vj) = (for i ^ j). The 

A: 



nr 



+ Amf/2 and 



mass eigenvalues corresponding to and za are, respectively, (m^) 2 
(m^) 2 = mf - Amf/2, where Amf ~ m^rriM <C Am^ latm . 

In the pseudo-Dirac scenario the probability of the v a — > up oscillation among the active 
neutrinos (a, /3 = e,n, r) can be written as: 



a/3 



E^ 



e j + e j 



(2.3) 



where J7 a i are the elements of conventional 3x3 neutrino mixing matrix (PMNS matrix). 
The phases for a neutrino with momentum k at the Earth is given by 



<l> 



± 



to 



1/2 



dt' , 



(2.4) 



where a is the time-dependent scale factor appearing in the metric (i.e., g^ u = (—1, a 2 ,a 2 ,a 2 )). 
The time of neutrino emission and the present time are respectively denoted by t and to. 
For baselines <C 100 Mpc, a(t)/a(to) — > 1 and we therefore obtain the standard result: 

$f = yjk 2 + (mf) 2 (t - t ). For the baselines 4irE u /Am 2 ol < L < 100 Mpc, the oscilla- 
tory terms originating from interference of different j (i.e., oscillation between pairs of states 
induced by Am^ and Am 2 tm scales) will average out and from Eqs. (2.3) and (2.4) we obtain 



P 



a/3 



(2.5) 



where Am 2 



(rrij) — (m- ) is the tiny mass-squared difference within the j-th pair of active 

Am? = Am 2 



i.e.. 



and sterile states. From now on, we assume equal Amj for all the pairs; „ 
for j = 1,2,3; we however keep the symbol Am 2 with its subscript "j". We will discuss 
the consequences of relaxing this assumption in the last section. By using Eq. (2.5) in the 
analyses of oscillation experiments, it is possible to extract or derive limits on the parameter 
Am 2 . Notice that the oscillation length corresponding to Am 2 is given by 



L, 



4vr£„ 



80 kpc 




(2.6) 
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which means that for the solar neutrinos with E v ~ MeV, oscillation length will be of the 
order of 1 a.u. for Amj ~ 10~ 12 eV 2 . The detailed analysis of the solar neutrinos gives the 
strongest bound on Am 2 , which is Am 2 < 1.8 x 10~ 12 eV 2 at 3a level [10]. For neutrinos 
coming from cosmological distances ~ Gpc with energy E v ~ 10 PeV, the baseline is much 
larger than the oscillation length (L 3> L osc ) for Am 2 > 10~ 16 eV 2 . However, for such 
distances the exact formulae in Eqs. (2.3) and (2.4) have to be used instead of Eq. (2.5). 
After averaging out the oscillatory terms induced by the solar and atmospheric mass-squared 
differences, the oscillation probability in Eq. (2.3) can be written as 

Pafi = W aj \ 2 \U Pj \ 2 cos 2 (^p-) where A$, = $t - *t . (2.7) 

j=l V / 

Notice that as long as we are in the ultrarelativistic limit with k 3> m f, the phase difference 
A$j remains proportional to Am 2 . Thus for Am 2 — > 0, we recover the standard formula 
for complete loss of coherence: P a p = Yli \U a i\ 2 \Upi\ 2 . From Eq. (2.4), the phases can be 
calculated as: 




(2. 



Thus, 



^^-t-^f-fM, (2 . 9) 

where E u = k is the neutrino energy at Earth which is related to the neutrino energy 
at source by E v = E®/(1 + z). Let us define a(io) = 1- Using a(t) = 1/(1 + z) and 
\dz/dt\ = H (l + z)y/n m (l + z ) 3 + n A , we obtain 

Am 2 f z dz' 

A$ 7 = -Dh / ; = , (2.10) 

3 2E U H J (l + z') 2 ^ m (l + z') 3 + ^A 

where Dh = c/Hq is the Hubble length with Ho = 71 km s" 1 Mpc" 1 denoting the present 
time Hubble expansion rate. Q m = 0.27 and £l\ = 0.73 are the matter and dark energy 
density of universe, respectively. Thus, for the phase difference appearing in Eq. (2.7), we 
obtain 

Am- 

= -^L{z) , (2.11) 

where 

r z dz 1 
L(z) = D H / = . (2.12) 

Jo (l + z')V^m.(l + ^') 3 +^A 

The red curve in Figure 1 shows the distance L(z) as a function of redshift. For com- 
parison, we have also shown the comoving distance d c (z) in blue [see Eq. (2.15)]. For a given 
E v , the asymptotic behavior of L(z) at high redshifts, sets a lower limit on Am 2 that can be 
probed by this method. This limit is of the order of 10 -17 eV 2 (£^/PeV). From this formula, 
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Figure 1: The red curve shows the distance L(z), defined in Eq. (2.12), as a function of 
redshift. The blue curve corresponds to the comoving distance defined in Eq. (2.15). The 
vertical axis is in the unit of Dh = c/Hq = 4.15 Gpc. 



we observe that for the purpose of probing smaller Am 2 , lower energy neutrinos have to be 
employed. A suitable choice can be supernova neutrinos with E v of the order of 10 MeV 
which can be sensitive to Amj down to 10 -25 eV 2 . Discussion about the neutrino oscillation 
in curved space-time can be found also in [11]. 

Let us now discuss two cases L L osc and L ~ L osc one by one. 

• In case that the baseline L is much larger than -L osc , the phases A<i>j will be very large. 
Since the spectrum of neutrinos from cosmological sources is continuous and the energy 
resolution of the detector is finite, regardless of the spatial distribution of the sources, 
the oscillatory terms induced by Am 2 will average out. Taking the average of the cosine 
term in Eq. (2.7), we obtain 

p^ = \Y,\ u ^\ 2 \ u ^\ 2 ■ ( 2 - 13 ) 

i 

Thus, for the case L 3> L osc we expect a suppression factor 1/2 of the neutrino flux 
with respect to the standard oscillation result due to the presence of tiny mass-squared 
differences Am 2 predicted in the pseudo-Dirac scenario. In particular ^2gP a p = 1/2, 
which reflects the fact that within the pseudo-Dirac scenario, half of the active neutrinos 
are converted to sterile neutrinos over large enough distances. 

In the limit L 2> L osc , this factor of half suppression turns out to be very robust. In 
particular, as we discuss below, possible matter effects at the source do not change 
it. If the effective potentials due to the matter effects (i.e., Vc = \f2Gpn e and Vn = 
— \/2Gi?n n /2) are much larger than Am? /2E U , the effective active-sterile mixing will 
be suppressed and the conversion from active to sterile neutrinos will be therefore 
suppressed. In order for the matter effects to play a significant role in the oscillation, 
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the following two conditions have to be fulfilled: (1) j Vn dt ~ J Vc dt ~ 1 or > 1; 
(2) Vn ~ ~ Am?/2E V (which is < Am^/2^). The first condition, which has 
been studied in detail in [12], is independent of the neutrino mixing and mass splitting. 
In a region where the second condition is fulfilled, there will not be any chance of 
conversion between vf and (with i ^ j) due to matter effects, as the splitting 
between m^~(~ m 4 ~) and m^" (~ mj) are given by (Am^ or A.m 2 atm ) / '(mj + m^) and 

are therefore too large. However, due to matter effects VLi — { v t + v 7)j\f2 can be 
converted to cos Qivf + sin 6ie ta ' za - , where #j and cxj depend on the matter profile in 
the medium through which the neutrino is passing. Suppose after exiting the matter, 
the neutrino traverses a large distance in vacuum. Then, 



P(v Li -)• v Li 



cosd^ + sm6 i e ia 'e i 



V2 



Thus, in the averaging limit P{i>Li v Li) = (cos 2 #j+sin 2 6i)/2 = 1/2. Now, consider a 
general state z/j n t = ^ SiU^i with £V s 2 = 1. The z/j nt state might correspond to a flavor 
state or might be a coherent mixture of them after passing denser regions. Suppose Vi n t 
arrives at such a region and transverses it and then travels a long distance L 3> L osc 

in vacuum. The oscillation probability to will then be P{yi n t v t-i) — Si I 1 2 if- 
This is the same results (suppression factor of two) that we would obtain without the 
matter effects on the active-sterile conversion. 

For the baselines comparable with the oscillation length, L ~ L osc , the cos 2 (3>j/2) in 
Eq. (2.7) does not necessarily average to 1/2 for an individual source and it can be 
smaller, leading to a strong suppression of active neutrino flux. However, for a flux 
of neutrinos originating from various sources with different baselines, one should also 
average over the different baselines. The suppression will then depend on the redshift 
distributions of the sources. In particular, if the neutrino sources are clustered around 
one or a few values of redshift z, for a given E u , the suppression factor can be made 
arbitrarily small by choosing the appropriate range of Am 2 . In order to quantify 
the significance of the suppression of neutrino flux, we define the effective suppression 
factor, S c g, as a function of Am 2 and neutrino energy at the detector, E u , in the 
following way: 

S cS (Am 3 ,E u ) i+ — , (2.14) 

where the sum is over the individual neutrino sources atredshifts z k . dN u /dE®(z k , E*>) 
is the energy spectrum of neutrinos from a source at z k and at the energy E® = 
E u (l + Zk). The symbol {X) show the average of the quantity X over the energy 
resolution of the detector. The factors of (1 + Zk) in the numerator and denominator 
come from the redshift of the energy interval: dE Q v = dE v {l + z k ). The factor l/d 2 c {z k ) 
is the comoving distance which accounts for the geometrical factor which enters in the 
calculation of neutrino flux: 

f z dz f 

d d z ) = D H / -==J—== ■ (2.15) 



y/n m (i + z>)3 + n A 
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For a fixed detected energy of neutrinos, Eq. (2.14) gives the suppression factor of 
neutrino flux, which depending on the redshift distribution of sources, can take values 
even less than 1/2 for an appropriate value of Amj. However, the neutrino telescopes 
such as IceCube, search for the neutrino flux from cosmological source over a range 
of energies. Thus, it is also convenient to define the overall suppression factor in the 
energy range (El,E%) as 



S eS (Amf,El,EZ) 



The quantity S e g, for a fixed value of Amj, gives a measure of the neutrino flux 
depletion in the whole energy range (E„,E„). 

3 Implications for Neutrino Flux from GRBs 

As discussed in the previous section, depending on the value of Am 2 and the baseline of the 
source(s), the flux of neutrinos from cosmological sources at the Earth can be depleted to 
different degrees. In this section we consider the specific case of GRBs as the sources of high 
energy neutrinos {E u ~ 100 TeV— 10 PeV). It is shown in [12] that inside the GRB sources 
J Vc,Ndt <C 1, so the matter effects inside the source are not relevant for neutrino oscillation. 
On the other hand, inside the Earth we shall have Vc,n 3> Am^/2E V , so again the matter 
effect will not be important because of the very tiny effective mixing angle. 

During the data-taking period of IceCube experiment (with two phases of 40 and 59 
strings out of the completed 86 strings implemented), about 30 GRBs have been observed 
by their gamma ray emission. From the X-ray observation of the afterglow it is possible 
to determine the redshift of the GRBs, which proves the cosmological origin of them [13]. 
However, this information is not available for all the GRBs. Figure 2 shows the redshift 
distribution of the observed GRBs for IC59. The red stars show the measured redshifts and 
blue crosses indicate the assumed values of redshift (by the IceCube collaboration) for GRBs 
without a measurement of redshift. Although for the majority of GRBs the redshift is not 
measured, obviously we can assume z > 1 for most of them; which implies the baselines 
L > Gpc. 

Following the discussion in section 2, for L S> L osc , the expected flux of neutrinos 
from GRBs will be reduced by half in the pseudo-Dirac scenario. For the neutrino energy 
E u ~ 1 PeV, this happens for Amj > 10 -16 eV 2 . Thus, taking into account the present upper 
limit from solar neutrinos [10], mass splitting in the range 10 -16 eV 2 < A?n 2 < 10~ 12 eV 2 
will deplete the GRB neutrino flux by a factor of 1/2. 

Now we discuss the interesting case of L ~ L osc . For the GRBs with redshift z ~ 1 and 
neutrino energies ~ PeV, this happens for 10" 18 eV 2 < Am 2 < 10~ 16 eV 2 . Following the 
discussion in section 2, we first calculate the suppression factor defined in Eq. (2.14) for GRB 
sources. However, as we depicted in Figure 2, for the majority of GRBs during the data-taking 
period of IceCube, the redshift is unknown. So, in order to calculate the suppression factor, 
we take two approaches: i) We calculate S e g just for the GRBs with a measured redshift. 
During the IceCube data-taking time, the redshift has been measured for 48 GRBs. ii) We 
assume that redshift distribution of GRB sources follows that of Star Formation Rate (SFR). 



COS" 



dN v {z k ,E°) (1+Zfc) , p 



I El 



dN v (z k ,E°) (l+z k ) 



~~dW 



\dc(z k )Y 



dE v 



(2.16) 
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Figure 2: Redshifts of the GRBs observed by their gamma ray emission during the data- 
taking period of IC59. The red stars show the measured redshifts and blue crosses indicate 
the assumed values of redshift by the IceCube collaboration. Data taken from [14]. 



Although in principle at least there is a delay in GRB rate with respect to SFR, to a good 
approximation the GRB redshift distribution is the same as distribution of core-collapsing 
massive stars which itself is proportional to SFR distribution (see [15] for the details of GRB 
and SFR redshift distribution correspondence). A good fit to the SFR data is given by [16] 



^) = C ^T^ — lYTW 2 — ' ( } 

where ip*(z) is the SFR per unit comoving volume and C is a constant (in the unit of 
Mq yr" 1 Mpc~~ 3 ). Thus, the number of GRBs with redshift (z, z + dz) is given by ip*(z) ■ dV c , 
where dV c is the comoving volume element given by ^itD nd 2 c (z)dz j \J £l m {\ + z) 3 + In 
order to calculate S c fi and S c s in Eqs. (2.14) and (2.16), we randomly generated 300 GRBs 
according to the distribution function tp*(z) ■ dV c . For the energy spectrum of GRB sources, 
we assumed that the sources at various redshifts have the same luminosity and same spectrum 
oc (E®)~ 2 ; i.e., we take dN u (z, E®) = N (E®)~ 2 dE® where N is independent of z. To obtain 
the average of cos 2 (A<£j/2) term in Eqs. (2.14) and (2.16), we assume a Gaussian distribution 
of energy with the width E u /W to model the energy resolution of the IceCube detector. 

Figure 3 shows the suppression factor S e g in Eq. (2.14) as a function of Amj for a fixed 
neutrino energy at Earth E u = 1 PeV. The blue curve corresponds to the SFR distributed 
GRBs and the red curve is for GRBs with measured redshift. For Am| < 10 -18 eV 2 , which 
corresponds to L < L osc even for the farthest GRBs, the suppression factor S c ff is equal 
to one, which means that the active neutrinos still do not oscillate to sterile states. For 
Arrij > 10~ 16 eV 2 , which corresponds to L ^ L osc even for the nearest GRBs, the suppression 
factor is equal to 1/2 which have been discussed in Eq. (2.13). For the intermediate values 
10~ 18 eV 2 < Am] < 10~ 16 eV 2 , the averaging discussed in Eq. (2.13) do not happen for all 
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1Q -18 1Q -17 1Q -16 1Q -15 

Am 2 (eV 2 ) 

Figure 3: The suppression factor S e f[ defined in Eq. (2.14) as a function of Am 2 . The 
blue curve corresponds to the SFR distributed GRBs and the red curve is for GRBs with 
measured redshift. The neutrino energy at Earth is fixed to E v = 1 PeV. 



the GRBs and thus the suppression factor can have values different from 1/2. Specifically, 
for Am 2 ~ 2 x 10~ 17 eV 2 , the flux of neutrinos from SFR distributed GRBs is completely 
depleted. The maximum depletion is slightly milder for the GRBs with known z which are 
more uniformly distributed. 

The position of dip in Figure 3 changes with E v . Thus, the strong suppression of the 
flux at a specific energy does not necessarily imply strong suppression in the whole energy 
range of observation in IceCube. We therefore calculate the overall suppression factor in 
Eq. (2.16) to estimate the depletion in the whole energy range. Figure 4 shows the S c q in 
Eq. (2.16) as a function of Am 2 for the energy range (E^,E^) = (0.1,3) PeV, corresponding 
to the energy range of IceCube experiment [4]. The color labels are the same as Figure 3. As 
can be seen, for Am 2 ~ 10~ 17 eV 2 , the overall suppression factor S e g can reach ~ 1/3, with 
slightly milder suppression for the sample of GRBs with measured redshift. 

In summary, by assuming appropriate value of mass splitting, the active-sterile oscil- 
lation in the pseudo-Dirac scenario can lead to strong suppression of GRB neutrino flux. 
The suppression can completely deplete the flux for some specific values of neutrino en- 
ergy, but in the whole energy window of the observation in IceCube it can reach to ~ 1/3. 
For the case Am 2 > 10~ 16 eV 2 , the energy spectrum of the GRB neutrino flux is not 
distorted and, independent of the neutrino energy, suppression factor is 1/2. But, for 
10~ 18 eV 2 < Am 2 < 10 -16 eV 2 , the S e s will have a strong dependence on neutrino energy. 
To illustrate this point, in Figure 5a and 5b we show the dependence of S c h on neutrino 
energy for two values of mass splitting Am 2 = 2 x 10~ 17 eV 2 and 5 x 10~ 17 eV 2 , respectively. 
In each plot, the blue and red curves respectively correspond to the SFR distributed GRBs 
and to the GRBs with measured redshift. As can be seen, by changing the value of Am 2 
the energy at which strong suppression takes place will change. Comparing the red and blue 
curves in each plot demonstrates that the precise position of the dip in the curves (corre- 
sponding to strong suppression in that energy) depends on the redshift distribution of the 
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Figure 4: The suppression factor S e s denned in Eq. (2.16) as a function of Am 2 . The 
blue curve corresponds to the SFR distributed GRBs and the red curve is for GRBs with 
measured redshift. The assumed energy range is (El,E%) = (0.1,3) PeV. 
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Figure 5: The suppression factor S e s defined in Eq. (2.14) as a function of neutrino energy 
For two values of mass splitting Am 2 = 2 x 10~ 17 eV 2 and 5 x 10~ 17 eV 2 . In each plot the 
blue curve corresponds to the SFR distributed GRBs and the red curve is for GRBs with 
measured redshift. 



GRBs. However, the presence of such a peculiar dip in the energy spectrum, independent of 
the precise position of it, can be interpreted as a hint for the presence of tiny mass splitting 
Am 2 . In principle, by reconstructing the energy spectrum of GRB neutrinos in the case of 
future observation and measuring the GRBs redshift distribution, it can be even possible to 
determine the values of Am 2 . 



4 Implications for Supernova type II Neutrinos 

Type II supernovae are the well-known and established sources for neutrinos beyond the solar 
system. The average energy of the supernova neutrino flux is ~ 10 MeV [17]. If a supernova 
explosion takes place within a distance of ~ 50 kpc, the neutrino detectors (including Super- 
Kamiokande [18], IceCube [19] and the proposed liquid Argon experiments [20]) can detect 
a few thousands of events in a few seconds (pointing towards a point source in case of the 
possibility to reconstruct the directionality). From Eq. (2.6), for supernova neutrinos the 
oscillation length will be L osc ~ 8 x 10~ 7 kpc(10 -12 eV 2 /Am 2 ) which is much smaller than 
the galaxy size or the distance between the only observed supernova up to now (SN1987A) 
and the Earth. Thus, even for Am 2 as small as 10 -18 eV 2 , we have L 3> L osc and the flux 
of neutrinos will uniformly be suppressed by a factor of two in the whole energy range of 
spectrum. As we discussed in section 2, in the averaging limit the suppression factor of 1/2 
will be independent of the possible matter effects inside the supernova. 

The expected neutrino flux from a supernova explosion depends on various parameters 
of the exotic medium inside the explosion region. Based on the supernova explosion models 
that are also verified by recent numerical simulations, the neutrino emission from supernova 
can be described by two phases of accretion and cooling. In the accretion phase the total 
flux of emitted neutrinos scales with MT acc , where M and T acc are the accreted mass and 
the temperature of medium, respectively. In the cooling phase the normalization of flux 
scales with i? 2 where R c is the neutrino-sphere radius. Considering the data available from 
SN1987A, the following values of these parameters and their la errors can be obtained [21]: 

M = 0.22±°;f| M , T acc = 2.4±g;2 MeV , R c = 16±§ km . 

It is easy to see that suppression factor of 1/2 for the SN1987A neutrino flux can be accom- 
modated within the error of normalization, so it will be a challenging task to investigate an 
energy independent suppression factor 1/2 in the total expected number of events. 

In principle, the same suppression happens also for the flux of relic supernovae neutrinos 
(or Diffuse Supernovae Neutrino Flux, DSNF). The DSNF comes from the contribution of 
all the supernovae explosion at various redshifts, during the history of universe after the era 
of star formation (see [22] for a review on DSNF). Up to now, no evidence for such diffuse 
flux is detected in the experiments. The strongest current upper limit on the flux comes from 
the Super-Kamiokande analysis of data collected over ~ 8 years [23], which is (at 90% C. L.) 
~ 1.2 cm" 2 s" 1 on the v e flux for E Ve > 19.3 MeV. Although the DSNF has not yet been 
observed, the sensitivity of the current experiments is very close to the theoretical predictions 
and the forthcoming neutrino detectors (such as Hyper-Kamiokande [24] and LENA [25]) or 
new methods of detection (such as doping the Super-Kamiokande with Gd [26]) would be 
able to detect this flux with a good significance. 

Theoretically, the DSNF can be calculated by integrating average neutrino flux from a 
typical supernova weighed with supernova occurrence rate over redshifts z = to z ~ 4 which 
corresponds to the beginning of star formation. In fact, the main contribution comes from 
the supernovae with z < 1. As discussed above, the normalization of a single supernova flux 
itself has a large uncertainty. On the other hand, the supernovae rate, which is calculated 
directly from the supernova observations or indirectly from the SFR, also has a relatively 
large uncertainty. It is shown in [27] that by taking into account these uncertainties, the 
theoretical prediction of DSNF can change by a factor of ~ 4 at 90% C.L. Thus, it can be 
argued that even the observation of DSNF in forthcoming experiment cannot rule out the 
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pseudo-Dirac scenario strictly. However, this scenario can "rule in" the models or supernova 
parameter ranges that predict a flux higher than the present or forthcoming experimental 
bounds. As discussed before we expect a distortion of the energy spectrum of DSNF for 
Am) ~ 10~ 25 eV 2 which increases the chance of identifying the active to sterile conversion 
despite the uncertainties in the total flux. Such small Am 2 will have no effect on the high 
energy neutrinos accompanying GRBs. 

5 Conclusions and Discussions 

Within the pseudo-Dirac scenario for neutrino masses, in addition to the two measured mass- 
squared differences (Am 2 ol and Am 2 tm ), there are tiny mass splittings in each pair of active 
and sterile neutrinos with Am 2 < 10 -12 eV 2 , where the upper limit comes from the analysis 
of solar neutrinos. In principle, these tiny values of mass-splitting can be probed by neutrinos 
coming from cosmological sources. The baseline of the neutrinos from cosmological sources 
with redshift z ~ 1 reaches ~ Gpc, which means that Amj ~ 10~ 17 eV 2 (E u /PeV) can be 
probed by the neutrinos with energy E v . However, due to the curved nature of space at large 
scales, the baseline saturates to its maximum value at z ~ 1 and for neutrinos from sources 
at higher redshifts, the baseline remains the same. Thus, the Am) ~ 10" 17 eV 2 (£„/PeV) 
can be interpreted as the smallest mass-splitting that can be probed by cosmological sources 
of neutrino with energy E v . 

We have shown that the active-sterile neutrino oscillation induced by Am) can partially 
explain the non-observation of the expected neutrino flux from GRB sources in the IceCube 
detector. The upper limit on the neutrino flux from GRBs derived by the partially deployed 
IceCube detector is ~ 3.7 times stronger than the expected flux, assuming that GRBs are 
the sources of UHECRs. For 10~ 16 eV 2 < Am) < 10 -12 eV 2 , the oscillatory terms induced 
by Am 2 will completely average out and the neutrino flux from GRBs will be suppressed by 
a factor of 1/2 , independent of the neutrino energy. For 10" 18 eV 2 < Am 2 < 1(T 16 eV 2 , 
the averaging does not take place for all the GRBs and the suppression factor can be smaller 
than 1/2 . For this case, we have shown that the average suppression factor over the whole 
energy range of the IceCube experiment can reach as small as 1/3, both when we assume 
that the GRBs are distributed according to the star formation redshift distribution and when 
we use the redshift distribution of GRBs for which the redshift is measured (which are less 
than 20% of the total number of GRBs studied during the IceCube data-taking period). 

With the present knowledge about the explosion mechanism in GRBs, we cannot con- 
clusively rule out /in the pseudo-Dirac scenario for 10" 16 eV 2 < Am) < 10~ 12 eV 2 . But, for 
10 -18 eV 2 < Am) < 10 -16 eV 2 , the suppression factor depends on the neutrino energy. As 
a result, the energy spectrum of GRB neutrinos will be distorted depending on the value of 
Am 2 . If the IceCube detector observes a flux of neutrinos associated to GRBs, by analyzing 
the shape of the reconstructed spectrum it will be therefore possible to test this scenario. 

For the whole range of 10 -18 eV 2 < Am) < 10~ 12 eV 2 , the oscillatory terms induced 
by Am) completely average out both for the point source and diffuse supernova neutrinos. 
Thus, the flux of supernova neutrinos will also be suppressed by a factor of 1/2, independent 
of the neutrino energy. Due to the large uncertainties in the theoretical prediction of flux nor- 
malization, investigating pseudo-Dirac scenario by supernova neutrinos will be challenging. 
For Am) ~ 10~ 25 eV 2 , we expect a distortion of the supernova neutrinos energy spectrum 
whose observation can be a more conclusive hint for the pseudo-Dirac scenario despite the 
uncertainties in the total flux. 
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Within the pseudo-Dirac scenario each active neutrino can be written as the sum of 
two quasi-degenerate mass eigenstates with equal contributions. However, if the scenario is 
extended such that each active state is the sum of n quasi-degenerate states, in the limit that 
all the oscillatory terms average out, we obtain a suppression factor of 1/n. 

Throughout the analysis in this paper, for simplicity, we took all the three Am 2 equal. 
In fact they can have different values. As long as they are all larger than 10~ 16 eV 2 , the 
oscillatory terms induced by each splitting will average out and the flux of all the active 
neutrino flavors will be suppressed by the same factor of 1/2. If some of Am 2 are smaller 
than 10 -18 eV 2 , the suppression of the neutrino flux from GRB sources will be reduced, as 
the corresponding state will remain active up to the Earth. If some of the Am 2 happen to 
be in the range 1CT 12 eV 2 < Am 2 < 1(T 16 eV 2 and different fr om each other, the distortion 
of the energy spectrum will be more complicated than the case that they are all equal. 
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